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Abstract: By using the AdS/CFT correspondence we study the deep inelastic scattering of 
an 7^-current off a = 4 supersymmetric Yang-Mihs (SYM) plasma at finite temperature and 
strong coupling. Within the supergravity approximation valid when the number of colors is large, 
we compute the structure functions by solving Maxwell equations in the space-time geometry of 
the AdS^ black hole. We find a rather sharp transition between a low energy regime where the 
scattering is weak and quasi-elastic, and a high-energy regime where the current is completely 
absorbed. The critical energy for this transition determines the plasma saturation momentum 
in terms of its temperature T and the Bjorken x variable: Qs = T/x. These results suggest 
a partonic picture for the plasma where all the partons have transverse momenta below the 
saturation momentum and occupation numbers of order one. 
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1. Introduction 

Over the recent years, there has been increasing evidence, coming from the experimental re- 
sults at RHIC and their theoretical interpretations [Q, ^, |^, |^ , and also from theoretical studies 
of the QCD thermodynamics [§, |6|, |^, that the hadronic matter produced after a high-energy 
heavy ion collision may interact rather strongly, in spite of being in the deconfined phase of 
QCD and having a relatively high partonic density. For instance, the success of theoretical ap- 
proaches based on hydrodynamics [^, |8|, which assumes local thermal equilibrium and vanishing, 
or small, viscosity, in describing collective phenomena like elliptic flow |^], suggests rapid 
thermalization and a low viscosity-to-entropy ratio for the matter produced at RHIC, which are 
hallmarks of a nearly-ideal fluid, with strong interactions. Also, the experimental results for 
the 'jet-quenching parameter' at RHIC |11, 12|, which is a measure of the rate at which highly 
energetic partons loose energy in the surrounding medium, have been interpreted lO, 113] to 



yield values which are too large to be explained by weak coupling calculations |15, [1^ (but this 
interpretation is not universally accepted; see, for instance, Furthermore, lattice studies of 

the QCD thermodynamics give evidence for a strong coupling behaviour (like the persistence of 



meson-like bound states ||18|, |19|, gO| ^ and strong deviations from the pressure of an ideal gas of 
quarks and gluons [§, ^) up to temperatures a few times the critical temperature for deconfine- 
ment. Such conclusions are corroborated by analytic calculations for the quark-gluon plasma 
showing that the weak-coupling expansion is too poorly convergent to be useful in practice for 



all temperatures of interest |22|, 23, 24 
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These and similar observations have urged the need for non-perturbative studies of the 
hadronic matter at relativistically high temperatures and densities. While lattice gauge theory 
is a privileged tool to non-perturbatively address static properties like the thermodynamics 
or the screening masses, its extension towards dynamical problems, like transport phenomena, 
dispersion relations, or the high-energy scattering, remains prohibitively complicated, and new 
methods are therefore required to systematically address such problems at strong coupling. The 
AdS/CFT correspondence (2^, although so far limited, in its most convincing formulation, to 
gauge theories which are 'simpler' (in the sense of having more symmetries) than QCD, is the 
most promising candidate in that sense. 

This method can most easily deal with the large-A'^ limit, with N the number of colors, 
where the gauge coupling g is small but the 't Hooft coupling X = g^N is large, in which case the 
N = 4: supersymmetric Yang-Mills (SYM) theory can be mapped onto a weakly-coupled string 
theory, that can be studied via semi-classical techniques. Leaving aside the structural differences 
between the = 4 SYM theory, which is conformal, and real QCD — these differences can be 
argued to be less important in well-chosen physical regimes, and, besides, some of them can 
be incorporated into extensions of the AA = 4 SYM theory (for which, however, the AdS/CFT 
correspondence is less firmly established) — , it is still not clear whether the aforementioned 
parametric conditions can be made consistent with the situation in QCD, where = 3 and 
(7 ~ 0(1) (giving A ~ 3H-6) in the interesting physical regimes. But even if a detailed, quantitative, 
comparison to real QCD (in particular, to the experimental data) would be premature, it is 
nevertheless clear that the AdS/CFT approach can provide valuable information about the 
non-perturbative behaviour of gauge theories, which should allow us to better constraint the 
physical reality of QCD from the strong-coupling end. 

Given these promising features, and the experimental imperatives at RHIC or LHC, it is not 
surprising that, over the last few years, there was a profusion of applications of the AdS/CFT 
techniques to problems of interest for high-density QCD. Following early applications to thermo- 
dynamics 1 27, and the pioneering calculation, by Policastro, Son, and Starinets, of the shear 
viscosity [29, 30 1, there was an intense activity towards computing the jet-quenching parameter 
1 31, 32, 33, 34 1, the energy-loss of a heavy quark |^, 36, 37, 38 1 or of a quark-antiquark pair 
|39, 4^, 41, 42, 43, 44], the diffusion rate for a heavy quark |45, the energy disturbances due 
to moving quarks |47, the Debye screening mass |49, the production rate for photons and 
dileptons |pl|, or the Bjorken expansion and the approach towards thermalization l53, 53, 54, pq] 



— all of that in the context of the strongly-coupled N = A SYM plasma at finite temperature 
(sometimes extended to include a chemical potential). 

Several of the studies mentioned above have been concerned with the long-range (Ax ^ 
1/T) or large-time (Ai ^ ^/T) behaviour of the strongly-coupled plasma, as relevant e.g. for 
hydrodynamics, thermalization, or transport phenomena. On the other hand, in order to study 
the propagation of 'hard' (i.e., highly energetic and relatively small) probes through the plasma, 
like jets or electromagnetic probes, it is essential to have a good understanding of the plasma 
structure on short space-time separations ^ 1/T, much alike the parton picture in perturbative 
QCD. Of course, at strong coupling there is a priori not clear whether the notion of a 'parton' — 
in the sense of a point-like constituent which behaves as quasi-free during the interaction with 
the external probe — makes sense in the first place, neither if such a 'parton', in case it exists. 



- 2 - 



should belong to an individual 'quasiparticle' (a thermal excitations with energies and momenta 
of order T), or rather is a property of the plasma as a whole. In other terms, is the partonic 
distribution of the plasma (again, assuming that this exists) the direct sum of the respective 
distributions for the constituent quasiparticles, with appropriate thermal weights, or rather is 
this qualitatively different ? 

Such questions are extremely difficult and below we shall not attempt to answer them in 
full generality. In particular, it is not yet understood whether a strongly-coupled gauge plasma 
admits a quasiparticle picture on the thermal scale 1/T, so like the Landau theory of a Fermi 
liquid, or the quasiparticle structure of the quark-gluon plasma emerging from resummations of 
perturbation theory |Q, ^. Fortunately, however, there is no need to properly understand the 
structure of the plasma on this scale 1/T so long as we are merely interested on the corresponding 
structure on much shorter space-time scales <C l/T. Indeed, the latter can be directly measured 
(at least, in a Gedankenexperiment) by an external probe with high energy and momentum 
{uj, q ^ T). From the experience with QCD we know that the most convenient measurement of 
that type — that whose results are most directly related to the parton structure of the target 
— is the 'deep inelastic scattering' (DIS) of a leptonic probe off the plasma. 

DIS at strong coupling in the context of the AdS/CFT correspondence has been so far 
considered |57, ^ only for the case where the target is a single hadron (a 'dilaton'). In this 
approach, the 'electromagnetic' probe which initiates the scattering is the conserved current 
associated with a particular U{1) symmetry (the '7?.-current'), whose associated '7^-charge' is 
also carried by the light degrees of freedom which are present inside the hadrons. By computing 
the current -current correlator in the hadron wavefunction, one can extract the same information 
about hadronic structure functions that would be obtained by DIS via a 'photon' coupled to 
the U{1) current. In terms of the standard kinematical variables and x, with = q^— uP' 
the virtuality of the current and x « Q^/s (at high energy s » Q^), the DIS structure function 
F2{x, Q'^) is a measure of the number of partons which carry a longitudinal momentum fraction 
X and occupy an area ~ l/Q^ in the transverse, impact parameter, space. 

In this paper, we shall use the same general setup — the scattering between the 7^-current 
and the plasma — to compute the structure functions of a strongly-coupled = 4 SYM plasma 
at finite temperature. It turns that, in this case, the formalism is quite different — in fact, 
somewhat simpler and also conceptually clearer — than in the case of a single-hadron target 
considered in Ref. |57, 58 1. There are several reasons for such differences: 

First, the string theory dual of the = 4 SYM plasma is unambiguously known, at the 
level of the original formulation of the AdS/CFT correspondence [26] : this is a 'black-hole' 
(more precisely, a non-extremal black three-brane; see Sect. 2 below for details) in a curved 
space-time geometry which is asymptotically AdS^ x S^. By contrast, in order to accommodate 
a hadronic state, the AA = 4 SYM theory (which has no confinement) must be 'deformed' in the 
infrared, in such a way to break down conformal symmetry. This deformation is not unique and, 
besides, its dual analog in the string theory is generally ambiguous. 

Second, the interplay between the large-A^ limit and the high-energy limit turns out to 
be much more subtle for a single-hadron target than for a plasma. This is in turn related to 
an essential feature of the strong-coupling problem, which is the deep connection between the 
distribution of partons and the issue of unitarity in DIS at high energy. As explained in Ref. 
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, at strong coupling, most of the partons are concentrated in the kinematical region where the 
scattering is strong and the unitarity corrections are important (the analog of the 'saturation', 
or 'color glass condensate', region of perturbative QCD [^). This is an important point that 
we shall try to motivate here via general arguments, and for which the subsequent calculations 
in this paper will provide an explicit realization: 

One can heuristically understand this point by extrapolating the picture of parton evolution 
in perturbation theory: Partons at large x tend to radiate and thus drop down at smaller values 
of X. At weak coupling, the emitted partons are predominantly soft (i.e., they carry only a tiny 
fraction x' ^ 1 of the longitudinal momentum of their parent partons), so, even for very high 
energies, there is still a substantial fraction of the parton distribution at relatively large values 
of X. These large-x partons carry almost all of the hadron energy and momentum, but they are 
unimportant for high-energy scattering, which is rather controlled by the bulk of the distribution 
at small-x. At strong coupling, on the other hand, there is no penalty for the hard emissions; the 
distribution of the energy among the child partons after a branching is essentially democratic, 
and hence the overall distribution can very fast degrade, via successive branchings, down to 
very small values of x. One therefore expects the structure functions at strong coupling to be 
concentrated at small values of x, but the question is, how small ? These functions are, of course, 
constrained by energy-momentum conservation — the small-x partons must carry the overall 
energy and momentum of the hadron — , but this constraint (a 'sum-rule' on F2) is not sufficient 
to determine the parton distribution. A more severe constraint comes from unitarity: in the 
kinematical region where the scattering is strong, in the sense that the scattering amplitude 
has reached the unitarity bound, the structure functions are, by definition, large, and hence the 
partons exist. Thus, at strong coupling, the search for the parton distribution is tantamount to 
understanding the unitarity problem for DIS. 

This is where the large~A^ limit becomes important: the elementary scattering amplitude 
is suppressed"*^ by a factor 1/N'^, so for a single-hadron target and in the strict large- limit 
{N ^ cxD at fixed energy), the scattering can never become strong, and thus the bulk of the 



partons cannot be seen. This is the situation considered in Ref. |57], and indeed it has been 
found there that the dilaton has no point-like constituents except at extremely small values of 
x (for a given resolution Q^), within a kinematical domain which squeezes exponentially to zero 
when increasing A. 

But the partonic structure of the dilaton reveals itself after relaxing the large- A limit, as 
we did in Ref. ||5^. Namely, we found that, for sufficiently large Q^, the partons are all located 
in the strong-scattering region at x < Xs{Q^), where Xs{Q'^) — h? / (N'^Q'^) is the 'saturation 
line' (a line in the kinematical plane (x,Q^) along which the elementary amplitude is constant 
and of order one) and A is the infrared cutoff which fixes the size of the dilaton. Moreover, 
the phase-space distribution of these partons turns out to be remarkably simple (and somehow 
reminiscent of the gluon distribution in the 'color glass condensate' at weak coupling |[59|]): 
there is essentially a single parton of a given color per unit cell in phase-space. This a posteriori 
legitimates the use of the 'electromagnetic' current as a probe of the parton distribution: in spite 
of the coupling being strong, the current can interact only with one parton at a time, and thus 



^There is no similar suppression for the DIS structure functions because the strength of the 7?.-current 
increases Hke A'^^, due to the color degrees of freedom of the fields which make up the current. 
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it can faithfully measure the parton number. Since, moreover, partons with very small x <^ Xg 
carry only little energy and momentum, it is clear that the hadron total energy and momentum 
is concentrated in the partons near the saturation line x = XgiQ"^)- 

Returning, after this long digression, to the plasma problem of current interest, we note 
that in this context one can simplify the problem by using the large-iV approximation without 
loosing the salient features: for a plasma target, the scattering can be strong even in the large- 
limit, because the plasma involves A^^ degrees of freedom per unit volume — as manifest from 



the fact that its entropy density scales like N^T^ 28] — , which compensates for the 
suppression of the elementary scattering amplitude. Note that, at this point, and at several other 
places in the paper, we use a heuristic language in which the plasma thermal degrees of freedom 
are treated as 'quasiparticles' with typical energies and momenta of order T, and the overall 
scattering process is viewed as the sum of elementary scatterings between these quasiparticles 
and the 7?.-current. This language, inspired by the situation at weak coupling, is admittedly 
ambiguous at strong coupling, and is used here only to gain more intuition into mathematical 
manipulations which by themselves are free of any ambiguity. 

Specifically, in the large-A^ limit of interest, the scattering between the 7^-current and the 
plasma can be described in the supergravity approximation, as the propagation of the gravita- 
tional perturbation induced by the current in the background metric of the black three-brane. 
The current-current correlator relevant to DIS is then computed from the action evaluated on 
the solution to the classical wave equation — an imaginary part in this solution being synony- 
mous of inelasticity in the scattering of the current. The wave dynamics is non-trivial in only 
one dimension — the radial dimension of AdS^, which plays the role of an 'impact parameter' 
between the current and the black hole. The relevant wave equation can be formally rewrit- 
ten as a Schrodinger equation in one spatial dimension (actually, two such equations, for the 
longitudinal and transverse waves, respectively). Then, the dynamics is controlled by the po- 
tential in this equation and, more precisely, by the competition between two important terms: 
a 'repulsive' term proportional to which by itself would keep the wave at the boundary of 
AdS^ (far away from the horizon of the black hole), and an 'attractive' term, proportional to 
the energy times the temperature, which tends to pull the wave towards the black hole. We thus 
distinguish between two physical regimes: 

(i) At relatively low energy and/or low temperature, such that x ^ T/Q, the repulsive 
term dominates, and the wave remains confined near the boundary. (For DIS off the plasma, 
X ~ Q'^/qT, and we recall that T/Q <C 1 for the physical problem of interest.) In this regime 
the scattering is weak and quasi-elastic (the imaginary part in the classical solution is extremely 
small, since generated via tunneling through the potential). Correspondingly, the DIS structure 
functions are exponentially small, e.g., F2 ~ exp{—c{xQ/T)^^'^}, a result that we interpret as 
the absence of point-like constituents in the SYM plasma having x ^ T/Q. 

(ii) At high enough energy, such that x < Xs{Q) — T/Q, the attractive term in the 
potential takes over, and then the wave escapes inside the bulk of AdS^, until it gets absorbed 
by the black hole. This absorption generates a large imaginary part in the solution, and hence a 
large contribution to the structure functions for DIS, which for x ~ is evaluated as F2 ~ N'^TQ 
(see Sect. 4 for more general results). These results for the structure functions represent the 
unitarity limit for the current-plasma scattering, which in this case is saturated by the complete 
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absorption of the current — a genuine 'black disk' limit. 

The physical interpretation of these results at small x in terms of partons in the plasma 
requires some care: the plasma being infinite, one needs to take into account the finite duration 
of the interaction, and also make a boost to a Lorentz frame where the notion of a parton makes 
sense (all the other calculations being done in the plasma rest frame). But after this is properly 
into account, it becomes clear that our results have a natural partonic interpretation, which 
moreover is consistent with the corresponding picture for a hadron, as obtained in Ref. [^ ]. 
Namely, for a given resolution Q^, the partons exist only at sufficiently small values of x, such 
that X < T/Q, and are homogeneously distributed in the three-dimensional phase-space, with 
occupation numbers of order one. Equivalently, for a given value of x, partons exist only at 
transverse momenta smaller than, or equal to, the saturation momentum Qs{x) — T/x. This 
value for the saturation momentum is consistent with the representation of the = 4 SYM 
plasma as an incoherent superposition of thermal quasiparticles. 

It is finally interesting to notice a similarity between our above estimate for the plasma 
saturation momentum and some results in the literature |4C, 41| for the screening length Ls{v, T) 
of a heavy quark-antiquark pair moving at velocity v in the hot = 4 SYM plasma. The 
screening length is the maximal separation for which the quark and the antiquark can be still 
connected by a string 'hanging down' in the radial direction of AdS^. In Refs. |pO| , |4l| ], one 
found Ls{v,T) ~ k(1 — f^)^/^/T with k a numerical constant (at least for v close to 1). Now, 



in the analogy with our DIS problem, the 'quark-antiquark pair' of Refs. |4y, 41] corresponds 
to the SYM system emerging from the 7^-current, which has a typical transverse extent 1/Q 
and a rapidity q/Q. It is therefore natural to identify our variables 1/Q and q/Q with the size 
L and the Lorentz gamma factor 7 = l/^/l — v'^ of the qq pair, respectively. In particular our 
saturation momentum Qs{x,T) should be compared to the inverse screening length l/Ls{v,T). 
To that aim, it is preferable to rewrite the result in Refs. |40, 41] as ~ 7^^; after replacing 
1/Ls Qs and 7 q/Qs, this translates into Q^{q, T) ~ which is parametrically the same 
as our result for Q<j, as alluded to above. It would be interesting to explore this correspondence 
in more detail. 



2. General setup and basic equations 

Following the general strategy with the problem of deep inelastic scattering, our objective will 
be to compute the retarded current-current commutator 

RA<1) = ^ j d^xe-^''-^^(xo)([J^(x),J,(0)]), (2.1) 

whose imaginary part determines the DIS structure functions. In the present context, the density 
J^{x) which enters Eq. ( |2.lD refers to an 7?.-current — the conserved current associated with 
a gauged U{1) subgroup of the 5'C/(4) global 7^-symmetry — , and the expectation value is 
understood as a thermal average, over the statistical ensemble corresponding to a = 4 SYM 
plasma at temperature T. The operator J^{x) for the 7?.-current receives contributions from the 
fermionic and scalar fields of the AA = 4 SYM theory. Accordingly, and following the example 
of perturbative QCD, we expect the imaginary part of Rfj,y{q) to give us information about the 
constituents of the finite-temperature plasma, just as the structure function of the proton gives 
information on its (partonic) structure in perturbative QCD. 
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In the limit where the Yang-Mihs coupling is small but the 't Hooft coupling \ = g^N 
is large, the AdS/CFT correspondence allows one to evaluate Eq. (2.1) in terms of classical 



supergravity in the metric of the AdS^ x black hole. The corresponding metric reads 

ds2 = ^ll^{-f{u)de + d^2) + -^—du" + R^d^l , (2.2) 

where T is the temperature of the black hole (the same as for the = 4 SYM plasma), R is 
the common radius of AdS^ and S^, t and x = {x,y,z) are the time and, respectively, spatial 
coordinates of the physical Minkowski world, u is the radial coordinate on AdS^, dOg is the 
angular measure on 5^ and /(u) = 1 - u^. Note that our radial coordinate has been rescaled 
in such a way to be dimensionless: in terms of the more standard, dimensionfull, coordinate r, 
it reads u = (^o/r)^, with rg = ttR'^T. Hence, in our conventions, the black hole horizon lies at 
u = 1 and the Minkowski boundary at u = 0. 

In order to evaluate Eq. (2.1), one needs to study the metric perturbation induced by the 



i?-current around the background metric (|2.2| ). The relevant gravitational wave is a vector 
field Afn{t,x,u) in AdS^, which obeys the classical equations of motion with given boundary 
conditions at u = 0. (Here, m = or u is the coordinate index on AdS^, with /i = 0, 1,2,3 
referring to a Minkowski coordinate.) Once the corresponding solution is known, the tensor i?^;^ 
can be extracted from the classical supergravity action evaluated as a functional of the boundary 
fields A^{t, X, 0) (see below). We shall assume the external current to be weak, so that the metric 
perturbations be small and the corresponding equations be linear in Am- Accordingly, we only 
need the supergravity action to quadratic order in Am, which reads (see, e.g., []64| ) 



647r2i? 



J d^x du V=5 ^'"^S"^ FmnFp, , (2.3) 



where Fmn = dmAn — dnAm, dm = d/dx"^ with = {t,x,u), and g = det{gmn)- The classical 
equations of motion generated by the action ( p. 3D are the Maxwell equations in the geometry of 
the AdS^ black hole. We shall work in the gauge Au = and choose the incoming perturbation 
as a plane wave propagating in the z direction: = (cj, 0, 0, q). Then we can write 

A^it,x,u) = e-'^'+"^'A^{u) (2.4) 

with the fields Afj_{u) obeying the following equations (below, i = 1, 2) 

■ujA'q + kfA'^ = (2.5) 



+ ^-jK + ^^A, = (2.6) 



A'o - ■^{k'^Ao + wkAs) = (2.7) 
uf 



where a prime on a field indicates a li-derivative and we have introduced dimensionless, energy 
and longitudinal momentum, variables, defined as 



Denoting a(u) = Aq{u), Eqs. ( |2.5[ ) and ( |2.7D can be combined to give 



a H H -p^ — a = 0, (2.9) 



which win be one of our key equations in what follows (the other one being Eq. ( |2.6D for Ai). 
The above equations (|2.5|)-(^^ have aheady been presented in the hterature (see, e.g., 



Refs. [ pO| , 61, 63, |6^), but in relation with other physical problems, corresponding to 
physical regimes very different from ours. These equations must be solved with the condition 
that the fields take generic values = An{u = 0) at the AdS^ boundary u = 0. Then Eq. (|2.7| ) 
implies the following boundary condition for a{u) 

limjna'('u)] = k{kAo + tuAs)]^^^ = k'^AiiO) . (2.10) 

For the solutions to be uniquely specified, an additional boundary condition is still needed. 
Following Refs. 64|, we shall require the solution to be a purely outgoing wave near the 



horizon at u = 1, where by 'outgoing' we mean a wave which is impinging into the black hole 
(and thus is departing from the Minkowski boundary). Physically, this corresponds to the fact 
that a wave cannot be reflected by the black hole, but only absorbed. Notice that, in the zero- 
temperature case where there is no black hole and the coordinate u extends to inflnity, the 
corresponding boundary condition is simply that the fields be regular at u ^ cxd. 

Once the classical solution is known, the next step is to compute the corresponding, 'on- 
shell', value of the action. Starting with Eq. ( p.3| ) and using the equations of motion to perform 
the integration over u, it is straightforward to deduce 



16 



J d^x [(Ao + J A;) a{u) - AduAiiu)] , (2.11 



where we have dropped a contribution coming from n = 1 in accordance with the prescription 
in Ref. |61, Note that the appearance of the factor in front of S is merely a consequence 



of our definition of the variable u (which scales like T^, so 9^ ~ 1/T^). If one returns to the 
dimensionfull radial coordinate r, then there is no apparent factor T^, and indeed Eq. ( |2.11| ) has 
a non-trivial limit as T ^ 0, corresponding to the vacuum polarization tensor for the 7^-current 
(see Sect. 3). Given the plane-wave structure in Eq. (0), the action density in Eq. (|2.1l|) 
is independent of = {t,x), so it is convenient to separate out the volume of space-time: 
S = J d'^xS. From the action density S, the tensor R^uiQ) is finally obtained as 

where A^ = A^{u = 0). Note that the ensuing tensor has mass dimension two, as it should. 

The tensor i?^;^ can be given the standard tensorial decomposition (see Appendix A), which 
shows that there are only two independent scalar components, Ri and R2, whose imaginary 
parts determine the two DIS structure functions, Fi and F2. (The precise definitions are given 
in Appendix A.) Since in practice we shall solve second-order differential equations with real 
coefficients, cf. Eqs. ( p.6| ) and ( |2.9[) , it is interesting to understand how an imaginary part in the 
respective solutions (and hence a non-vanishing contribution to the DIS structure functions) 
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can arise in the first place. This is generated via the aforementioned boundary condition near 
u = 1, which allows for the absorption of the gravitational wave by the black hole. 

For a given temperature T of the target plasma, the scalar functions Ri, or Fi, with i = 1,2, 
depend in general upon two kinematical invariants, that we shall conveniently choose as the 
virtuality of the T^-current and the Bjorken x variable. These are defined as 

q2 =q2_^2 ^ ^ ^ (2.13) 

^ ^ ' -2{q ■ n)T 2ujT ' ^ ' 

where is the four-velocity of the plasma in a generic frame, and the second expression for x 
holds in the plasma rest frame, for which n'^ = (1,0,0,0). Unless otherwise specified, in what 
follows we shall always work in the plasma rest frame. We shall consider the large-Q^ and 
high-energy kinematics, where ^ ^ and hence uj q. These conditions allow for 
both small (x ^ 1) and large (x ~ C(l)) values of x, but in what follows we shall be mostly 
interested in small-x regime where q ^ Q"^ /T. 

To conclude this section, let us present an alternative form for our key equations, ( |2.6| ) and 
(|2.9| ), which is more insightful and also better suited for constructing approximate solutions via 
WKB techniques. Via simple manipulations, these equations can be brought into the form of 
the (time-independent) Schrodinger equation in one spatial dimension, that is, ip" — Vip = 0. 

Consider first Eq. ( [2.91 ): when rewritten for the new field ip{u) = \/u{l — u'^) a{u), this takes 
the form (with K"^ = k"^ — zu'^) 



1 



A(i + 6n2-3n4)+i^2_^2^2 

4-0, 



^^ = 0, (2.14) 



which is of the Schrodinger type, as anticipated. In the interesting regime at > > 1, the 
potential V{u) in ( 2.14| ) is well approximated by 



V 



u(l — v?)"^ 



1 +i^2_^2^2 



4u 



(2.15) 



This describes a potential barrier, whose shape is illustrated in Fig.l, and also in Fig. 2a, for 
three different physical situations, corresponding to different regimes for the ratio k/K^ : (i) 
k/K^ < 8/(3V^) in Fig.la, (ii) k/K^ = 8/(3^3) in Fig.lb, and (iii) k/K^ > 8/(3^3) in 
Fig.lc. As it should be clear from Fig.lb, the critical value k/K^ = 8/(3\/3) corresponds to the 
case where the height of the potential vanishes at its peak. Note that a value of 0(1) for the 
ratio k/ corresponds to a value x ~ T /Q <C 1 for the Bjorken variable. 

We can understand much about the solution to ( p. 14 ) by inspection of these figures: When 



k/K^ < 8/(3\/3) (the situation at intermediate energies), there is a high potential barrier (cf. 
Fig.la), with classical turning points ui ~ l/(4i^^) and U2 — K/k. (Note that ^ n2 ^ 1 in 
the interesting regime where k ^ K ^ 1, with k <^ though.) We then expect the solution 
'4>{u) to be concentrated within the classically allowed region at n < Moreover, the DIS 

structure functions are expected to be extremely small in this case, since an imaginary part in 
the classical solution can develop only via tunneling through the high potential barrier. 

On the other hand when k/K^ > 8/(3\/3) (the high -energy case, cf. Fig.lc), there is no 
potential barrier any longer, so the gravitational wave can easily fiow into the black hole and 
thus get absorbed by the latter. We then expect a large imaginary part to R^i,. 
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(a) 



(b) 



(c) 



Figure 1: 

The potential V{u) in Eq. ( p^ ) for three values of the ratio k/K^: (a) k/K^ < 8/(3^3), (b) 
k/K^ = 8/(3^3), and (c) k/K^ > 8/(3^3). For the fi gures to look more suggestive, all the 
chosen values for k/K^ are relatively close to the critical value 8/(3\/3)- 




(a) 




(b) 



Figure 2: 

The potentials V{u) corresponding to longitudinal waves, cf. Eq. (|2.15|) (left) and, respectively, 
transverse waves, cf. Eq. ( 2.171 ) (right), are represented for several values of the ratio k/K^, 
corresponding to physical regimes well separated from each other. 



Similar conclusions apply to the transverse modes Ai as well, although the corresponding 
argument is slightly more involved, and perhaps less intuitive. With the substitution (j){u) = 
y^(l — u^) Ai{u) (for either i = 1 or i = 2, the respective equations being identical), Eq. ( |2.6| ) 
takes the Schrodinger-like form 

^" - n ^ 2N2 [^^ - ^^"^ -u]cP = (2.16) 

u[\ — w^y 

where for the present purposes the potential can be approximated by 

^ = ^T^-^ - ^'"1 • (2.17) 
ml — w^Y 
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This potential is illustrated in Fig. 2b, for two values of the ratio k/K^ . As manifest in these 
figures, the potential barrier is now concentrated near the boundary at ti = 0, within a distance 
u < I/K"^, whereas the classically allowed region (i.e., the region where V{u) < 0) starts at 
u = K/k. With increasing energy, the barrier does not disappear anymore, rather it gets 
squeezed towards u = 0, in such a way that its effects become smaller and smaller. At low 
energy, the wave can penetrate into the bulk only up to a small distance u ~ l/K"^ away from 
the boundary. But when the energy is so high that k/K^ ~ C?(l)i the penetration distance 
~ I/K"^ becomes of the same order as the classical turning point at K/k, and then the wave 
can freely escape in the allowed region at n > K/k, and thus get absorbed by the black hole. 

These simple observations will be confirmed and substantiated by the subsequent analysis 
in this paper. 



3. Low energy: the multiple scattering series 

In this section we shall consider the low-energy regime at k/K^ <^ 1, cf. Fig. la, where the 
effects of the term proportional to k'^ in the potential V{u) (in either Eq. ( |2.15| ), or ( 2. IT]) ) can 



be treated in perturbation theory. Note that, in terms of our original variables uj and q, cf. 
Eq. (p.8|), the condition k <C K^ amounts to qT'^ <C Q^- Hence, for a fixed virtuality Q^, the 
'low-energy' regime can be also understood as a low-temperature one, T <C {Q^ /q)^^'^ , and the 
perturbative expansion that we shall shortly construct can be alternatively viewed as a multiple 
scattering series, or a low-temperature expansion. 

Clearly, even when k <^ K^, this perturbative expansion cannot work for arbitrary values of 
u : when u > K/k, the energy-enhanced term cx k'^ in the potential becomes the dominant term 
there, which is responsible for the existence of the classically allowed region at u > K/k. Thus, 
not surprisingly, the perturbative treatment of the finite-energy/temperature effects cannot 
account for the contributions due to tunneling, which are genuinely non-perturbative and will be 
estimated in Appendix B within the WKB approximation. But if one leaves these contributions 
aside (they are exponentially suppressed anyway; see Appendix B), then perturbation theory 
should work reasonably well in the small-u region at u < 1 /K'^ , which is the relevant region for 
computing the 7^-current correlator, cf. Eqs. ( p.ll[ )-( p.l2D . 

The main result that we shall arrive at in this section could be characterized as negative: 
we shall find that for k <^ K^ the DIS structure functions are strictly zero when computed to 
all orders in the multiple scattering (or 'twist') expansion. But the subsequent analysis is still 
interesting in that it provides the twist expansion for the real part of Rfj,u- In particular, from 
the leading term in this expansion (the single scattering approximation), we shall be able to 
deduce a couple of energy-momentum sum rules which will be very useful later on. 

In the interesting region at u < l/K"^ <C 1, our key equations (l2.6|) and (^) simplify to 



and, respectively. 



^2 

A'' Ai = -k^uA, , (3.1) 

u 



1 K"^ 

a" + -a' a = -k'^ua . (3.2) 

u u 
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In writing these equations, we have separated the terms oc k'^ in the r.h.s., anticipating that 
they are going to be treated as 'smah perturbations'. For consistency with the present approx- 
imations, which ignore the phenomenon of tunnehng, the above equations must be solved with 
the condition that the fields vanish as u — > oo. (This would be the correct boundary condition 
in the zero-temperature limit T — > 0, and it remains the appropriate boundary condition for a 
perturbative treatment of the finite-temperature effects.) 

Consider first Eq. (|3.2| ); after a change of variable C = 2Ky/u, this becomes 

5p + ^5^-l)«(C) = -j^a(C). (3.3) 

The zero-temperature limit^ of this equation, that is, 

d2 1 d 



describes the (longitudinal) metric perturbation induced by the 7?.-current in AdS^ in the absence 
of the black hole (the supergravity dual of an 7^-current propagating through the gauge theory 
vacuum). The general solution to ( p.4| ) is a linear combination of the modified Bessel functions 
Kq and Iq. The coefficient of Iq is set to zero by the condition of regularity as C — *■ c«, while 
that of Kg is fixed by the boundary condition at C = 0, cf. Eq. ( p.lO| ). One thus finds 

aW(C) = -2A;2^z.(0)Ko(C). (3.5) 
The general equation ( p. 2] ) can be given a formal solution via Green's function techniques : 

poo /_^2a/4\ 

a(C) = aW(C)+y^ dC'GiCC') ^-j^j a{0, (3-6) 
with the Green's function G{(, (') obeying 

+ J^-i)g(c,C') = 5(C-0, (3.7) 

together with the following boundary conditions 

G(C,C') ^0 as C^oo, 
C^G(C,C') - as (3.8) 

It is easily checked that the corresponding solution reads 

G(c, C') = -C {Ko(C)io(C')©(C - CO + Ko(C')io(C)e(c' - 0} • (3.9) 

The 'solution' ( tj.6| ) is truly an integral equation for a(C), which generates the multiple scattering 
series through iterations — here, for the longitudinal wave. 



^The limit T ^ of the present equations may look tricky since we have defined the dimensionless variables 



in Eq. (2.8) by dividing though T. However, in the zero-temperature case, one can view T in Eq. (2.J) as an 



arbitrary reference scale, introduced in order to define dimensionless variables. This scale cancels out in the final 



results for the current correlator at T = 0, as one can check on the examples of Eqs. (3.17) and (3.1^) below 
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Consider similarly the transverse sector. By replacing At = Ai(0)(h{C), with = 2Ky/u, 
within Eq. (3.1), one finds 



dC^ C dC CV ^ ^ 

The zero-temperature version of this equation is solved by h^^^ = Ki(C), which obeys Ch^'^HC) ~^ 
1 as C — *■ 0, as it should. (The other solution Ii(C) is rejected by the condition of regularity at 
infinity.) The general equation ( p. 10 ) can then be rewritten as an integral equation similar to 



Eq. (3^) with a^^'^ — > h^'^^ and the following Green's function 

G(C,C') = -C'{Ki(C)ii(O0(C-O + ii(C)Ki(O0(C'-C)}. (3.11) 

As a simple application of the previous results, let us now compute the first two terms 



in the low-temperature expansion of the current -current correlator (2A) — that is, its zero- 



temperature piece which represents the vacuum polarization tensor of the 7?.-current, and 
the first temperature-dependent contribution R^f^u , which describes the scattering between the 
7^-current and the M = 4 SYM plasma in the single-scattering, or 'leading twist', approxima- 
tion. To that aim, we need the first two iterations in the above integral equations for a{u) and 
Ai{u), evaluated near u = (cf. Eqs. (pID-(^l2D). 

To the order of interest, we can write a{u) = a^^\u)+a^^\u), where a^'^\u) has a logarithmic 
singularity as u ^ 0, as expected according to Eq. (|2.10| ), 



,(0) 



k^AL{0){lnK'^ + \nu + 2-f + Oiu)) (3.12) 



(7 = 0.577... is Euler's constant), while a^^\0) is finite and equal to 

a«(0) = Al{0)-^ J dCC^K2(C) = --^Al{0). (3.13) 

As for the transverse fields Ai{u), these are needed up to linear order in u, i.e., to quadratic 
order in ( ; one finds Ai{u) = Af'\u) + a\^\u), with 

Af^\u) ^ Ai{0)\^l + uK^{lnK^ -l + lnu + 2-f)^ (3.14) 
A^\u) ^ Am^,. (3.15) 



When the T = fields in (|3T^ ) and (|1|) are used to evaluate the vacuum action S^'^\ cf. 



Eq. ( |2.11] ), the result exhibits a logarithmic divergence coming from the limit u — > 0. As usual 
in the AdS/CFT context, this singularity is interpreted as a ultraviolet divergence in the dual 
gauge theory, to be removed via renormalization. To that aim, it is important to return to the 
original variables r, u;, q, and = — uj"^ , to make it clear that the UV 'counterterms' are 
indeed temperature-independent. Then, the relevant terms in the action are 

\nK +\nu = In^^+ln = In^+ln^-^, (3.16) 
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where the T-dependence has disappeared, as anticipated, and A plays the role of the substraction 
scale on the gauge theory side. For convenience, we renormalize by dropping the last term in 
the above equation together with the finite term 27. We thus obtain S^^^ = J d^xS^^\ with 

= In 5^ [(^^0 + u^A.f - Q'At ■ At ] , (3.17) 

where we have introduced the transverse vector notation At = (j4i, A2)- From this expression, 
one can immediately deduce the vacuum polarization tensor (with r/^^^ = (—1, 1, 1, 1)) : 

R^^Hl) = {v,. - ^) R?\Q') with <(Q2) = ^ In . (3.18) 

This is transverse, as required by current conservation, and moreover it has exactly the same ex- 
pression as in zeroth-order (one loop) perturbation theory. This 'non-renormalization' property 
of the 7^-current polarization tensor has been already observed in the literature, and proven to 



be a consequence of supersymmetry [35|. 



Similarly, by using the finite-T contributions to the fields, Eqs. ( |3.13| ) and (3.15), one can 



compute the respective contribution to the on-shell action, 5^^-* = J d'^xS^^\ which is manifestly 
ultraviolet-finite : 



5(1) 



2^2t^4 „2 r 



30 



(3.19) 



u=0 



From (|3l9| ) one can determine the single-scattering, or low-temperature, part of the tensor 



This has the structure exhibited in Eq. (AT) with the following, 'leading-twist', expressions for 
the scalar components Ri and R2 ■ 

which are both real : as anticipated at the beginning of the section, the DIS structure functions 
vanish in the leading-twist approximation, and in fact to all orders in the twist expansion — 
indeed, all the terms generated by iterating the integral equation (|3.6| ) for the longitudinal field, 
or the corresponding equation for the transverse fields, are obviously real. 

Note the l/x"^ behavior of R^\ which is the hallmark of the graviton exchange and reflects 



the fact that the contributions to ii^j/ computed in Eq. ( 3.20 ) come from the twist-two and spin- 
two operator T^j^ (the energy-momentum tensor) in the operator product expansion (OPE) of 
the current -current correlator ( |2.1[ ). This is the only leading-twist operator which survives in 
the OPE at strong coupling, since the other twist-two operators with spins j > 2 acquire large 
anomalous dimensions ~ A^^^ 00. It is quite remarkable that the OPE coefficients oiT^^ that 
we have (indirectly) computed at strong coupling are exactly the same as the corresponding 
coefficients at weak coupling, as we shall demonstrate via an explicit zeroth-order calculation 
in Appendix C. This non-renormalization is yet another manifestation of the high degree of 
symmetry of the AA = 4 SYM theory (see also Ref . Q ) . 

Similarly, the multiple scattering series previously discussed can be interpreted as the ex- 
change of arbitrarily many gravitons. One simple way of understanding the lack of an imaginary 
part in these multiple graviton exchanges is to note that the gravitons carry no four-dimensional 
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space-time momentum, as reflected in the fact that the metric only depends upon the radial 
variable u in AdS^. Hence, because of energy-momentum conservation, the graviton exchanges 
cannot create on-shell final states, which would be the source for inelasticity. 

We thus conclude that in the intermediate energy/low temperature regime at /c <C K"^ 
(or, equivalently, at relatively large values x ^ T/Q for the Bjorken variable), the only non- 
trivial contributions to the DIS structure functions Fi oc Im Ri arise via tunneling through the 
potential, and thus are necessarily small. In Appendix B, these contributions will be estimated 
in the WKB approximation as Fi ~ exp{— c(i^^//!;)^/^}, where the prefactor c is a number of 
0(1). This estimate confirms that Imi?j remains extremely small so long as k <^ K^. We 
thus draw the rather striking conclusion that the strongly-coupled plasma has essentially no 
point-like constituents at x larger than Xs ~ T/Q. 

Finally, let us mention an interesting consequence of the leading-twist results in Eq. ( |3.20| ) , 
which will be very useful in what follows. Introducing the variable z = 1/x and assuming 
standard analytic properties for the current-current correlator in the complex z plane, one can 
relate the behaviour of Ri{z) near 2 = 0, where Eq. ( 3.20| ) applies, to the integral of the DIS 



structure function F^ oc Imi?j along the cuts on the real axis in the physical region at \z\ > 1. 
One thus obtain the following sum-rules (see Appendix A for details) 

£ = 18T^[ dxF2{x,Q^), (3.21) 
Jo 

£ = 45T^[ dxFLix,Q^), (3.22) 
Jo 



where Fi and F2 are defined in ( [A.2| )-( |Al3 ), Fl = F2 — 2xFi is the longitudinal structure 
function, and 

£ = ^ (3.23) 

is the energy density of the = 4 SYM plasma in the strong coupling limit: £ = 0oo, with 



Qp^y the energy-momentum tensor of the plasma, cf. Eq. (A. 4). The appearance of the energy 



density in the l.h.s.'s of equations ( 3.21j ) and ( |3.22 ) is in fact natural: as we shall further explain 



in Sect. 5, the integrals in their r.h.s.'s are proportional to the energy density carried by the 
plasma constituents, as probed in DIS with a resolution scale Q^; this should be the same as the 
total energy density in the plasma, and in particular be independent of — which is precisely 
the content of Eqs. (|3.21[ )-( |3l22|) . 



But the previous results in this section also show that the relatively large values of x, 
such that X > T/Q, give only tiny contributions to the structure functions, which die away 
exponentially at large and hence cannot ensure the fulfillment of the sum-rules. Therefore, 
the only way for these sum-rules to be satisfied is that the integrals in their r.h.s.'s be saturated 
by contributions from 'partons' at smaller values of x < T/Q. This corresponds to the 'high- 
energy' situation in Fig. Ic, to the analysis of which we now turn. 

4. High energy: deep inelastic scattering 

In this section, we shall consider the high-energy (k > K^), or small-x (x < T/Q), regime, 
which is the most interesting regime for our present analysis, since this is where the deeply 
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inelastic scattering truly occurs. In this regime, the potential barrier becomes ineffective — it 
has either completely disappeared (in the longitudinal sector, cf. Figs. Ic or 2a), or become so 
narrow that it gives no significant attenuation (in the transverse sector, cf. Fig. 2b) — , and 
then the gravitational waves induced by the 7^-current can propagate towards large values of 
u ~ C(l), until they reach the black hole horizon at u = 1 and thus get absorbed. As explained 
in Sect. 2, this absorption manifests itself via imaginary parts in the classical solutions, that we 
shall first compute, and from which we shall then deduce the DIS structure functions. 



By lack of exact solutions to the wave equations (2^) and (|2.9| ), we shall consider approxi- 
mations which are valid for very high energies, such that k ^ K^, but which cannot capture the 
transition from quasi-elastic to deeply-inelastic scattering, which takes place around k ~ K^. In 
Appendix D, we shall construct approximate solutions valid for generic values of u, by perform- 
ing piecewise approximations (in particular, the WKB approximation) and then matching the 
intermediate solutions with each other. Here, however, we shall use a simpler strategy to calcu- 
late the classical action (|2.11|). To appreciate this strategy, let us first recall what was the main 



difficulty with this calculation: although the action involves the classical solution near u = 



alone, as manifest on Eq. (2.11), this solution is generally sensitive to the dynamics at large 
u ~ 0{1), via the 'outgoing wave' boundary condition that one has to impose near the horizon. 
The important simplification that appears at high energy is that this boundary condition can 
now be imposed already at relatively small values u <^ 1, where the general equations reduce to 
simpler ones, that can be solved exactly. Indeed, in the absence of any potential barrier, there 
is no mechanism to generate reflected waves at intermediate values of n < 1; hence, an incoming 
wave cannot be tolerated in the solution not even at u <C 1, since it would necessarily describe 
reflection off the black hole. This argument will be confirmed by the more general construction 
in Appendix D, which will provide the same small-n solutions as obtained below in this section. 

Note an additional, important, simplification which occurs at high energy: when k ^ 
K^, the term in the potential involving the virtuality of the current becomes negligible as 
compared to the other terms there, for all the relevant values of u (for both longitudinal and 
transverse modes). This means, in particular, that our subsequent discussion also applies to a 
time-like {K"^ = k"^ — < 0) current, provided its energy is high enough {k ^ |-f^|^)- 

Indeed, consider the longitudinal sector first. When increasing u from n = 0, the last term 
oc /c^n^ in the potential ( p.l5| ) becomes comparable to the first term oc l/Au already at the very 
small value uq = l/{Ak'^)^/^ , at which the term oc is still negligible. Hence, the latter is never 
relevant, as anticipated. In particular, for n ^ 1, the potential simplifies to 



V ~ — — 



for u <C1 , (4.1) 



which has a peak at u ~ mQ) cf. Fig.lc. By performing the corresponding approximations on 



Eq. (|2.9|) , this equation becomes 

1 



a" + -a + k^ua = , (4.2) 



u 



which can be easily solved: after changing variable according to ^ = |/cm^/^, we obtain 



^ + ^^ + l)««)=0 (4.3) 
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which has the general solution 



«(e) = ciJo(0+C2No(0, (4.4) 

where Jq and Ng are the usual Bessel and Neumann functions. Recalling the behaviour of these 
functions near ^ = 0, one sees that the boundary condition ( p.lO| ) fixes the coefficient C2, 

C2 = ^^l(0), (4.5) 

but has no consequence for ci. To also determine the latter, we shall impose, as announced, 
the outgoing-wave boundary condition at sufficiently large values of u. Note that, although 
u is small, n <C 1, the argument ^ of the Bessel functions becomes large, ^ 1, for all the 
values of u far beyond the peak of the potential: u » mq ~ l//c^/^. In that region, one can 
use the asymptotic expressions for the Bessel functions, that is, Jo('?) — cos(.^ — 7r/4) 

and No(0 ^ sin(C - vr/4). If one also remembers the exponential factor yielding the 

time-dependence, cf. Eq. (|2.4|) , it becomes clear that a purely outgoing-wave solution a{t,^) oc 
Q-i{u}t-^) -g Qbtained by choosing 

ci = -ic2 . (4.6) 
One thus obtains the following expression for the longitudinal solution in the region u <C 1 

where H^^^ = Jo + i^o is a Hankel function. 

A similar discussion applies to the transverse waves, which satisfy Eq. (U), or ( ^) . In 
the high-energy regime at A; > K^, one can neglect the effects of the extremely narrow potential 
barrier located at < u < K/k. Indeed, the width K/k of the barrier is much smaller then the 
distance ~ l/K"^ over which the solution near u = would start to significantly differ from its 
boundary value at u = 0. In the small-u region at K/k < u ^ 1, the potential ( p. 17 ) reduces 
to y ~ —k'^u and then both Eq. ( |2.6| ) and Eq. ( p.l6| ) reduce to 

A'l + k'^uAi = . (4.8) 

This is an Airy equation whose general solution can be written as a linear combination of 
Ai(— li/cs) and Bi(— m/cs) or, equivalently [^], in terms of the Bessel functions of argument 
v = 1/3. We choose this latter representation, for more symmetry with the previous discussion; 
we thus write (with ^ = as before) 

Aio = [ci Ji/3(e) + c2Ni/3(o] , (4.9) 

where C2 is determined from the value of at u = and we again choose ci = —ic2, in order for 
the solution to become a purely outgoing wave when u ^ l/fc^/^. One finally gets the following 
result at small u : 

M^) ^ MO) (0 H« (^fc^'/') for n « 1 , (4.10) 
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which now features the Hankel function H^^g = J1/3 + 

By putting together the previous results ( ^^ ) and ( 4.10| ), one can evaluate the on-shell 
action according to (2.11); this gives S = J d^xS with 



48 



1 ^ 
7 + ln- 



ITT 



+ 



r2(i/3) 



(4.11) 



Note the emergence of the imaginary part in the action ( [4.11| ) , which has the right sign (Im S > 
0) to describe dissipation, i.e., to yield positive contributions to the DIS structure functions. 
A simple calculation using Eqs. (IID, O, (^) and ( [A3| ), finally leads to the following 
expressions for the structure functions at small x <^ Xs ^ T/Q : 



Fl 



16r2(l/3) V3 ) ' 
N^Q^x 



F2 - 2xFi 



967r2 



(4.12) 
(4.13) 



which represent our main result in this paper. Although strictly valid only for x <C x^, these 
results remain parametrically correct also in the transition region at x ~ x^. For x ^ x^, on the 
other hand, the structure functions are negligibly small, as discussed in Sect. 3. 

To render the above results more transparent, it is convenient to rewrite them in terms of 
the conventional variables for DIS, x and Q^, and to also introduce the transverse structure 
function Ft = 2xFi, such that F2 = Ft + Fl. Then Eqs. ( 4.12 ) and ( 4.13| ) imply the following 
parametric estimates: 

2.^^2T^ W^^/^ 



Ft{x,Q 



X \ T2 

r2 /x2Q2 



X 



2^2 



(4.14) 



which show that, in the very small-x regime at x <C T/Q, the longitudinal structure function 
is negligible as compared to the transverse one, Fl <C Ft, and thus, somehow surprisingly, 
an analog of the Callan-Gross relation applies: F2 — 2xFi. This looks surprising since it is 
quite different from what happens in the case where the target is a single hadron, at either 

where in the high-energy limit Fl and Ft are 



weak coupling'^ |59|, or strong coupling 
parametrically of the same order. 
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5. Saturation and the partonic structure of the plasma 

The results of the last two sections are conveniently described using Fig. 3 where r = Inl/x is 
the 'rapidity' and p = ln{Q^ /T"^). For a given p ^ 1 and values of r below the saturation line 

''in QCD at weak coupling, the Callan-Gross relation holds only in the Bjorken scaling regime at relatively 
large x, where the structure functions are dominated by the valence quarks and depend very weakly upon Q^. 
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r =ln- 




Quasielastic 



P = WlT 



Figure 3: Proposed 'phase diagram' for DIS off a JV = 4 SYM plasma at high energy and strong 
coupling. 



Ts{p) = p/2, meaning x Xs — T/Q, the structure functions are extremely small (cf. Eq. ( [B.lOl ) 
in Appendix B), 

Fl^t ~ N^Q'^x exp { - c{x/xs)^^^} cx exp { - ce^""^-^)/^} for r < Ts{p) , (5.1) 

while for values of r significantly above that line (x <C Xg) the structure functions take on the 
values given in ( [4.12 ) and ( 4.13| ). The transition between these two regimes when crossing the 



saturation line is expected to occur within a rapidity interval At ~ 0{1). 

The saturation line can equivalently rewritten as Ps{t) = 2t, and then the estimates (|5^ 
for the structure functions at small r are tantamount to 

Fi ~ exp { - c(Q/g,)i/2} ~ exp { - ce^"-"^^/^} for p > p,(r) , (5.2) 

with the saturation momentum 



Qlir) = T^e^^ = T^e^" = ^. (5.3) 

-|2 ^ /-)2 



Such a small value for Fi for large Q ^ QHt) is qualitatively consistent with previous calcula- 



tions of the dilaton structure functions at strong coupling 58 1 , although some quantitative 
differences remain. In these previous works, one has found that the higher-twist terms domi- 
nate the dilaton structure functions at large Q^, thus yielding a fast decrease with Q^, which 
is however power-like, Fi{x,Q'^) cx (l/Q^)^ with A > 1, rather than exponential as predicted 



by Eq. (5^). Some of these higher-twist contributions are naturally absent from the present 
analysis, since suppressed in the large-A^ limit. (This is the case for the diffractive processes 
considered in Ref. ISq], which in the present framework would correspond to multiple scattering 
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off a same ttiermal quasiparticle. Tlie corresponding scattering amplitude starts at order 
and tience it is suppressed at large N even after multiplication by the number ~ N"^ of thermal 
degrees of freedom.) On the other hand, the higher-twist contributions due to protected oper- 
ators, as discussed for a dilaton target in Ref. [57|, would survive in large-A^ limit, but they 
are removed by the requirement of energy-momentum conservation. (Being homogeneous in the 
four physical dimensions, the plasma cannot transmit any energy or momentum via a single scat- 
tering.) We interpret this smallness of Fi at relatively large x to mean that for x ^ Xg = T/Q 
there are hardly any point-like excitations (partons) in the SYM plasma. 

In what follows, we shall rather focus on the more interesting situation at x < x^, or large 
rapidity r > Ts{p), where partons do exist, as we shall see. A natural place to look for a 
partonic interpretation is at the level of the sum rules ( |3.21| ) and (3.22). By inspection of our 
previous estimates (4.14) for the structure functions, it is easy to check that (i) the integrals 
in Eqs. ( 3.21| ) and ( p. 22 ) are dominated by values of x of order Xs (for which the transverse 
and longitudinal structure functions are of the same order of magnitude), and (ii) the results 
of these integrations are of the right order of mag nitude, namely of 0{N'^T^), to ensure the 
fulfillment of the sum rules. For instance, for Eq. (p.21|) we can write 



£ = 18T^ [ dxF2{x,Q^) ^T^xF2{x,Q^) , (5.4) 

Jo x=T/Q 

where xF2{x, Q^) ~ N'^T'^ when x ~ T/Q, as manifest from Eq. ( [4. 141) (recall that F2 = Ft+Fl). 
Ours results in Eqs. (P^-(^l3|) are not accurate enough to also check the numerical coefficients 
in front of the sum rules (this would require a more precise study of the transition region at 
X ~ Xs). But the parametric estimates in Eq. ( 4.14 ) are sufficient for our present purpose, which 
is to develop a partonic picture for the strongly coupled plasma. 

Before we proceed, let us first recall the interpretation of the structure function F2 in the 
more familiar context of perturbative QCD. In that case, F2{x,Q^) is a dimensionless quantity 
interpreted as the quark distribution in the proton target, i.e., the number of quarks which are 
localized in impact parameter space within an area ~ 1 /Q"^ fixed by the resolution of the virtual 
photon, and which are distributed in longitudinal phase-space within a unit of rapidity (Ar ~ 1) 
around the rapidity r = ln(l/a;) fixed by the Bjorken variable. In what follows, we shall boldly 
propose a similar interpretation for the strongly-coupled plasma, and then critically examine 
the most sensible points in our proposal. 

Unlike the proton, or dilaton, structure functions, which are dimensionless, the plasma 
structure functions Fi as computed in this paper have dimensions of (area)~^. This makes it 
natural to try and relate these functions to the density of partons per unit area in the transverse 
plane (x, y) (the impact parameter space). The 7?.-current with » probes an area ~ l/Q^ 
much smaller then the typical area ~ 1/T^ covered by a 'thermal quasiparticle' in the plasma 
— i.e., a typical thermal excitation with energy and momentum of order T (in the plasma 
rest frame). This means that the current can see 'inside' a quasiparticle, and thus probe its 
elementary constituents, or 'partons'. More precisely, it will simultaneously scrutinize inside 
all the quasiparticles located within one coherence length in the longitudinal direction z. The 
notion of coherence length is particularly important for what follows, and so is also the choice 
of an appropriate Lorentz frame in which the parton interpretation makes sense. So, let us open 
a parenthesis at this point, in order to better explain these concepts : 
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(i) The partonic picture makes sense in a frame where the current has low energy and 
a relatively simple internal structure, so that it can act as a probe of the target. Here, it 
will be convenient to use the Breit frame where the 7^-current is a standing wave. Namely, 
if one boosts the plasma by an amount rj where coshrj = q/Q, then in this boosted frame the 
current has time and z-momentum components a;' = and q' = Q. This current is naturally 
absorbed by partonic constituents of the boosted plasma having momenta of order Q. Indeed, 
the partons participating in the collision have a longitudinal momentum fraction x, and thus a 
typical momentum ~x(Tcoshr/) ~(5 in the boosted frame. 

(ii) Before boosting, the current correlator ( p.l| ) is sensitive to longitudinal distances Az < 
2q/Q'^, as is suggested by writing the space-time dependence in the integral there as 

^-iu]t+iqz ^ ^-iq{t-z)+iQ'^t/2q ^ (5.5) 

where we have used ~ 2g(g — io) for q^ ^ Q^. This shows that the integration in ( |2.1| ) can 
extend over the coherence time Ate ~ '^q/Q'^, corresponding to a coherence length Azc ~ '^q/Q^ 
in the plasma rest frame. After the boost, this length gets Lorentz contracted (note that the 
current is decelerated) down to a value Az'^ ~ (2g/Q^)/ cosh?] ~ 1/Q. 

Let us now return to Eq. ( |5.4| ) and try to interpret this sum rule in the Breit frame. In 
the l.h.s., we would like to construct the energy density per unit area, dE' /d'^b, of the region of 
the plasma which is explored by the 7^-current in this boosted frame. As a component of the 
second-rank tensor Q^u, the (three-dimensional) energy density £ transforms in the boost by a 
factor (coshr/)^. The 7^-current probes a slice of the plasma with longitudinal extent Az'^ in the 
boosted frame. Hence, dE'/d'^b ~ {£ cosh^ rj) A z'^ ~ £{2q/Q'^) cosIit]. Multiplying both sides of 



(|5.4|) by {2q/Q ) cosh rj, one gets 



dE' „ , /l^, ^2 



- xTcoshi](-F2{x,Q^)] . (5.6) 

Jx=T/Q 

As before mentioned, the quantity xTcoshrj ~ Q in the r.h.s. is the longitudinal momentum of 
the constituent (parton) which interacts with the 7^-current. It is therefore natural to interpret 



1 



X 



dn 



x=T/Q d'^b 



-F2{x,Q') ^_~:j27 (5.7) 



x=T/Q 



as the number of partons per unit area within a longitudinal slice of the plasma, with the width 
of the slice equal to the coherence length (which is 1/Q in the Breit frame, and q/Q'^ in the 
plasma rest frame). This interpretation, which here has been inferred from the sum rule ( ^.4] ), is 
in fact natural in view of the standard partonic interpretation of F2 at weak coupling, as alluded 
to before. The only new feature with respect to the case where the target is a single proton^ is 
the factor 1/x in the l.h.s.: this is a Lorentz-invariant measure of the amount of matter in the 
plasma in the longitudinal slice explored by the current. Namely, this has been generated as 
(say, in the plasma rest frame): TAzc ^ 1/x, where T is the density of quasi-particles per unit 
length and Azc is the longitudinal extent of the interaction region. 



*There is no such a factor in the case of a single-hadron target since there the whole longitudinal extent of 
the hadron lies within one coherence length for the virtual photon. 
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For what fohows, it is useful to notice that the parton density in the r.h.s. can be equivalently 
written 

dn , dn dn dn 

d26 " diM26 - p'^dz'd% " d7d26 ' 

where Az' ~ 1/Q is the longitudinal extent of the slice in the boosted frame and p'^ ^ Q (the 
z-momentum of a struck parton) is the same as 1/Az', as it should by the uncertainty principle. 
In writing the last equality, we have identified the rapidity interval dr = p'zdz' . 

By using Eq. (5.7) together with the previous estimate (4.14|) for F2 = Ft + F^, one finds 



N^Q^ for x^T/Q, (5.9) 



drd26 



which, remarkably, has the same parametric form as in a weakly-coupled gauge theory , and 
hence it admits a similar physical interpretation. Namely, when interpreted in the Breit frame, 
Eq. ( |5.9D gives the total number of partons (per unit area) having a longitudinal momentum 
fraction equal to x (with x <T /Q) and with transverse momenta p_L < Q. Since this number 
appears to be of order N'^Q'^, we conclude that there is a number of order one of partons (of a 
given color) per unit of phase-space: 

1 dn 

(The factor dr plays the role of p'^dz', cf. Eq. (|5.8|) , so the above phase-space density is an 
occupation number, in the proper, three-dimensional, sense.) This is similar to pQCD in the 
sense that the parton occupation number saturates at sufficiently low transverse momenta, below 
a critical scale Qsix) which grows like a power of 1/x. In QCD, saturation is a reflection of 
unitarity in a corresponding scattering process. Where is the unitarity limit here? Viewed on 
the gravity side of the AdS / CFT correspondence the gravitational wave induced by the TZ- 
current is completely absorbed at the horizon of the black hole and that absorption takes place 
over a time less than or equal to the coherence time, 1/xT, of the wave. This is, in effect, a 
unitarity limit for scattering of the gravitational wave (very much similar to the corresponding 
limit for dipole scattering in the familiar dipole factorization for DIS at high energy ||59|] ). 

Recently, the phenomenon of parton saturation in relation with the unitarity limit for DIS 
has also been identified at strong coupling, in the case where the target is a single 'hadron' (a 
dilaton) [^]. Interestingly, our above result (^]^) for the saturation momentum of the plasma 



is consistent with the corresponding result for a single hadron in Ref. |58], once the assumed 
structure of the SYM plasma in terms of quasiparticles is taken into account. Namely, the 
quantity is proportional to the density of partons per unit area in impact parameter space. 
In the case of a single dilaton, Ref. [^] has found 

Qs{x) = — TTTT (one dilaton target) , (5-11) 

with 1/A a measure of the dilaton transverse size. When moving to the plasma, the dilaton 
gets replaced by thermal quasiparticles with individual size ~ 1/T. To account for the degrees 
of freedom relevant to DIS, one must sum over color (this yields a factor N'^) and also over the 
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number of quasiparticles within a longitudinal slice of width Azc ~ ^/Q^ in the plasma rest 
frame (which gives an additional factor TAzc ~ After replacing A — > T in Eq. ( |5.11| ) 

and putting these various factors together, we end up with the previous result, Eq. ( p.3\ ), as 
anticipated. This is an important check of the internal consistency of our proposed partonic 
description — it comforts our idea that DIS off the plasma at high ^> should measure 
the internal constituents of the thermal quasiparticles composing the plasma. This check is 
particularly non-trivial in view of the fact that the unitarization mechanisms at work appear to 
be very different in the two cases — disappearance of the potential barrier for the plasma case, 
respectively, diffractive scattering via multiple graviton exchanges in the dilaton case. 

We now turn to the case x <^ T/Q, which turns out to be quite subtle. Previously, we 
argued that the typical interaction time is of the order of the coherence time Ate ~ q/Q^ of 
the incoming current. This argument, however, ceases to be valid at very high energy, where 
the gravitational wave gets absorbed (reaches the horizon) on a time scale shorter than Ate- A 
heuristic way to understand this is to recall that, when the energy is so high that the barrier 
has disappeared, cf. Fig. Ic (namely, for qT'^/Q^ ^ 1), the -Fi structure function becomes 



independent of Q^, as manifest on Eq. ( 4.12 ). On the other hand, the definition ( p.l| ) of the 
current -current correlator involves an explicit dependence upon Q^, via the exponential factor 
inside the integrand, conveniently written as in Eq. { f).5\ ). The only way for this dependence to 
disappear at high energy is that the integral over t in Eq. (2.1) be cut at some time which is 



considerably shorter than the coherence time Ate- This requires the lifetime of the gravitational 
wave (before being absorbed by the black hole) to be shorter than Ate- For a given energy 
q, this lifetime can be estimated as q/Ql{q), since = Ql{q) is the smallest value of Q for 
which Eqs. (l2.9|) and (|2.6|) still have a Q-dependence. Here, Ql{q) is the saturation momentum 
expressed as a function of q, and is obtained from the condition qT/Ql = Qs/T (the condition 
to lie on the saturation line in Fig. 3) as Ql{q) = (gT^)^/^. Note that q/Ql{q) is indeed much 
smaller than q/Q"^ in this high energy regime. 

Because of this short lifetime of the high-energy gravitational wave, we believe that the 
quantity Fi{x,Q'^) ~ (l/2x)F2(x, Q^) (which, we recall, is the dominant structure function 
when X <C T/Q) is actually being determined by interactions of the 7^-current with partons of 
size 1/Ql{q), rather than with partons of size 

1/Q2. Indeed, Eq. (|1|) imphes 



Fi{x,Q^) = Fi{xsMl) where Ql{q) = {qT^f^ and Xs{q) = T /Qs{q) . (5.12) 

Let us give another argument leading to the same conclusion. We recall that in DIS in the 
QCD dipole picture, and in the rest frame of the target, the size of the dipole emerging from 
the electromagnetic current expands with time as Ax±_ ~ ^JTj^q |7^], so that it takes a time 
t ~ 2q/Q'^ for this dipole to reach a size Aa;_|_ ~ 1/Q- Assume that a similar estimate applies 
for the SYM system emerging from the 7?.-current (the analog of the QCD color dipole); then, 
after a time tg ~ 2g/Qg(g), which is the lifetime of this SYM system before being absorbed by 
the plasma, its size gets only as big as Ax_l ~ 1/Qs{q)-, which indicates once again that the 
partons at scale Qs{q) are the relevant degrees of freedom. 

We shall conclude this discussion, and also the paper, with a critical analysis of the main 
assumptions that we have made in reaching Eqs. (^]^) and ( ^.9]) — the equations at the basis 
of our partonic interpretation, (i) We have assumed the plasma to be made of constituents 
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('quasiparticles') having momenta on the order of T (in the plasma rest system) when measured 
on a resolution scale T. The fact that entropy density and energy density scale as N'^T'^ and 
N'^T'^ suggest that this is the case, but this understanding is, perhaps, not completely clear. 
Note that, for the present purposes, we did not need to specify the actual nature of these 
'quasiparticles', which at strong coupling would be a most difficult task, (ii) We have also 
assumed that the T^-current directly measures individual constituents at scale Q. In QCD the 
electromagnetic current provides such a measurement in the leading order renormalization group 
formalism. At next-to-leading order, ambiguities occur in separating the measured partons 
from the probe; however, these ambiguities are effects of order a(Q'^) and cannot affect general 
conclusions as to numbers of partons in a hadron or plasma. In SYM, we have taken the coupling 
large so that the separation between the probe and the partons to be measured is not sharp 
anymore. In reaching (5.7) and (|5.9|) we have assumed that, up to factors of order one, the TZ- 



current couples to individual constituents of the plasma and that this coupling is not strongly 
renormalized. Because of these subtleties we feel that our results have to be taken with caution, 
and that a deeper understanding of the partonic structure of the plasma in strong coupling SYM 
is highly desirable. 
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A. Structure functions: Definitions and sum rules 

In this appendix we remind the reader of the tensor structure of R^^ and derive a sum rule 
relating the expectation of energy momentum tensors in the plasma to deep inelastic scattering 



on the plasma. R^^, is defined in Eq. ( |2.1| ), The tensor structure must be given in terms of 
and the plasma four-velocity n'^ = {nt,nx,ny,nz), with n'^ = (1,0,0,0) corresponding to the 
plasma at rest. Then current conservation plus the symmetry property R^y{q) = Rufi{—Q) imply 
the following general structure (with tj^i, = ( — 1,1,1,1)) : 



R2 ■ (A.l) 



nf^Ui, - ^^{rif^qu + n^q^) + ^^lyj • ^lY 

The two scalar functions -Ri and R2 depend upon the two invariants and x introduced in 
Eq. (|1|), and they are even functions of x. We define the DIS structure functions as (note 
that n ■ q = —lo is negative) 

Fi = — luiRi, (A.2) 



27r 
-(n • q) 



2ttT 

Writing the energy-momentum tensor of the plasma as 

3 



Imiia • (A.3) 



= [VpLy + '^n^n^) - with 8 = , (A.4) 
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we can rewrite the leading-twist results in Eq. ( |3.20| ) as 

This rewriting makes it clear that the calculation of the leading-twist contribution to Rfj^i, in 
Sect. 3 amounts to computing the coefficients of the energy-momentum tensor in the operator 
product expansion for the current-current correlator. 

To deduce the sum rules ( 3.21 ) and ( p.22| ), we shall write z = 1/x and assume the standard 



analytic structure for the functions Ri{z) in the complex z plane. Namely, Ri{z) is an analytic 
function everywhere in the complex plane except for two cuts along the real axis (from z = —oo 
to 2 = — 1 and, respectively, from z = \ to z = oo). Then our previous results in Eq. ( |A.5| ) 
express the dominant behaviour of Ri{z) near z = Q. Using this information together with 



Eq. ([A. 51) , we can successively write 



4£: f dz R2 _^ f°° dz 2i Imi?2 



9Q^ J 27ri z Ji 27ri z 

= - [\x^-^ = ^ [' dxF,{x, Q') (A.6) 



TT Jq X 







where the contour in the first integral is a small circle surrounding the origin. This is then 
distorted in the complex plane in such a to wrap around the two branch cuts which give equal 
contributions. (We assume the integrand to vanish sufficiently fast as \z\ — > 00 to be able to 
neglect the contributions of the large circles closing the contour.) Similarly, 

dz Ri 2 dz 2ilmRi 



15r^ / 2Tri z^ Ji 27ri z'^ 

[ dxxlmRi=4 [ dxxFi(x,Q2) (A.7) 
Jo Jo 



2 

vr 



B. Low energy: tunnel effect 

In this Appendix we shall use WKB techniques to estimate the probability for inelastic scattering 
via tunnel effect in the intermediate energy regime at Q <C § <C Q^/T'^, where the potential 
barrier is high. The argument turns out to be non-trivial because the imaginary part of the 
classical solution — which, we recall, is the measure of inelasticity in the scattering — gets 
built via a 'double-tunnel effect' (see below), for which the WKB approximation is generally 
not reliable. Yet, as we shall later argue, in the present setup this approximation should be 
reliable for the imaginary part of the solution. 

We shall focus on the longitudinal wave (the corresponding discussion of the transverse wave 



is entirely similar) and use the wave equation in Schrodinger form, cf. Eqs. ( 2.14 )-( 2.15 ). We 
shall construct our global approximation for ■0 by matching approximate solutions valid in three 
different domains: (i) u close to zero, (ii) u inside the potential barrier in Fig. la, and (iii) 
relatively large u, on the right side of the barrier. As usual, the imaginary part in the solution 
will be generated by the condition that ip{u) be a purely outgoing wave at large u ~ 0{1). 
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(i) For relatively small u, the potential in Eq. (2.15) can be approximated as 



y ~ i 

u 



4u 



for 



< n <C K/k . 



Then the general solution can be written as: 



(B.l) 



(B.2) 



where the coefficient Ci is fixed by the boundary condition at u = 0, Eq. ( 2.10| ), as Ci = 
— 2A;^^l(0). This approximation is similar to the zeroth order perturbative solution ( p.5| ) in 
Sect. 3 except that, now, the coefficient C2 in front of Iq is allowed to be non-zero because of 
the different behaviour assumed at large u. The imaginary part, ImC2, of this coefficient is the 
quantity that we are primarily interested in, because this quantity determines, via Eq. ( |2.11| ), 
the imaginary part of the 'on-shell' action. 

(ii) For values of u inside the potential barrier, ui < u < U2 with ui ~ l/(4i^^) and 
U2 — K/k the two classical turning points in Fig. la, the solution can be constructed via the 
WKB approximation, which yields 

1 



C3 exp |- r du'^/n^l + C4 exp I r du'^/V{ 



(B.3) 



In applications of the WKB technique to the tunnel effect, the analog of the second term in the 
equation above is generally omitted, since beyond the accuracy of this approximation. However, 
in so far as the imaginary part of the solution is concerned — which, we recall, is our main 
interest here — the inclusion of this term is both essential and justified, as we shall later argue. 

The approximate solutions ( [B.2D and ( [B.3[ ) have a common validity range at ui < n ^ K/k, 
and thus can be matched with each other in this window. By also using the asymptotic behaviour 
of the modified Bessel functions, as valid for u ^ ui, one finds 



/vr 



C2 = 2^/7? C4 , 



(B.4) 



(iii) For u ^ U2, the WKB solution is similar to the one constructed in Sect. 4 and reads 

C 



tp{u) ~ 



exp 



i / duyJ-V{u') 



(B.5) 



where we have selected only the outgoing wave, i.e., the one propagating towards the black 
hole. (This corresponds to choosing C5 = in Eq. (|D.5[) .) The above coefficient C is the 
same as C4 in Eq. ( D.5D , but its precise value is irrelevant here (it would merely determine the 
normalization of the wave near u = 1, cf. Eqs. ( p.2| ) and (|D.7D ). Rather, what matters is 
the relative normalization of the coefficients C3 and C4 in the solution (|B.3| ) inside the barrier, 
which in turn is fixed by matching Eqs. (B.3) and (B.5) near u = U2- This matching cannot 
be done by directly comparing these two solutions, as they have no overlap with each other. 
Yet, the proper matching procedure is standard in the WKB literature |71, ^ (this requires a 
study of the exact behaviour near n2, which can be done by linearizing the potential and then 
recognizing the Airy equation), and here we shall simply list the result: 

C -j7r/4 ^ i 

2 



C3 



D 



DCe 



-i7r/4 



DC. 



(B.6) 
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where D is the WKB attenuation factor (in the usual context of quantum mechanics, this 
describes the decrease in the intensity of the wavefunction after passing the potential) : 



"2 



Ui 



D = e^p\-2 / du^/V{u) } . (B.7) 



By comparing Eqs. (B.4) and (|B.6|) one finds 



C2 = i^DCi, (B.8) 



which implies the following behaviour for (|B.2| ) near the boundary at u = (recall that ip{u) ~ 
y/ua{u) for small u) : 

a{u) ~ A;^^l(0) [ In + (In u + 27) - inD ] . (B.9) 

This is our main result here. It shows that, in the presence of a high potential barrier, the 
imaginary part of the solution near n = gets built via a double-tunnel effect. This is 'double' 
since the relative strength of the imaginary part versus the real part is D, and not ^/D. This 
result is in fact natural: this imaginary part is the feedback of the absorption taking place 
near u = 1 on the gravitational perturbation at the Minkowski boundary. First, the incoming 
perturbation, which is purely real, has to cross the barrier to approach the black hole, then, 
after the scattering takes place near u = 1, the imaginary part thus generated in the solution 
must propagate backwards and cross the barrier once again, before being measured (in the form 
of DIS structure functions) at u = 0. 

By using Eq. (|B.9| ) together with the corresponding equation for the transverse sector, one 
can finally compute the DIS structure functions generated through tunneling in this low-energy 
(or low-temperature) regime. One thus finds quite similar expressions for the longitudinal 
(Fl = F2 — 2xFi) and transverse {Ft = 2xFi) structure functions : 

i^. {i=L,T). (B.IO) 



It is easy to check that the integral in Eq. ( |B.7] ) is dominated by the region in u where the 
potential can be simplified to V{u) ~ {K^ — k'^v?)/u, which is the same as the potential ( |2.17| ) 
in the transverse sector and for u -^1. Then, the attenuation factor is essentially the same (to 
leading exponential accuracy) for both the longitudinal and the transverse waves, and can be 
estimated as 

D^e^v{-c{K^/kf'^] with 2r^(V4) (g-^^) 



This explains the estimates (5.1)-(^r2|) for Ft and Ft 



Let us finally explain why, in the present context, we think that it was justified to keep 
the second term in the WKB solution ( |B.3| ). Generally, this term is discarded in applications 



of the WKB method 72] since it is exponentially suppressed as compared to the first term 
there (recall that C4 oc -DC3, cf. ( |B.6D ), and hence it is much smaller than the corrections to the 
prefactor in that first term, which are only power-suppressed (in this case, by rational powers 



of k/K^). However, in the present problem, the first exponential in (B^) is matched onto the 
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real part, oc Ci Kq, of the solution at small u; hence this large exponential term is strictly real, 
and so would be all the higher-order terms, neglected by the WKB approximation, which would 
correct its prefactor. Accordingly, the second exponential in ( |B.3| ) is the only one which can 
develop an imaginary part, and this imaginary part is therefore correct to WKB accuracy. To 
conclude, the WKB approximation cannot be trusted for the real part of the coefficient C4 in 
(|B.3| ), but only for its imaginary part, which is the quantity of interest for us here. 

C. The operator product expansion at weak coupling 

In this Appendix we shall show that, when computed to lowest-order in perturbation theory, 
the coefficient of the energy-momentum tensor in the operator product expansion (OPE) of the 
current -current correlator (^]^) is exactly the same as the corresponding coefficient in the strong- 
coupling limit, as implicitly computed in ( ^.20D (or in Eq. ( |A.5| )). This nonrenormalization 
property reflects the high degree of supersymmetry of = 4 SYM (see, e.g., [^]). In the 
perturbative calculation of the OPE to follow, we shall keep only the operators which mix with 
the energy— momentum tensor T^^. 

In AA = 4 SYM, there are six scalars (pm in the vector representation of 50 (6) and four Weyl 
fermions ipi in the fundamental representation of SU{A). The IZ symmetry current corresponding 
to the generator = diag(l/2, — 1/2, 0, 0) is 

= ^(V^ia'^V'' - V^^aV) + \ {(peD^c/^^ - <P5D^<P6 + c/^iDf'cps - (psD^^'Pi) , (C.l) 

where = (1, — <t), <t being the Pauli matrices. By contracting fields with free propagators, it 
is straightforward to derive the relevant part of the OPE: 



^ m— zi 



^ \i=l,2 m=3,4,5,6 j 



i_ (^r^f^^rj'^P - ^^^/^ - + {l^i + l^'p + (™nglet terms)) 



1 / ^ 



) + (nonsinglet terms)) + • • • . (C.2) 



where in writing the second equality we have projected onto the S'C/(4) singlet operators and 
denoted 

6 6 

= E = E ^miDaiDpcl,^ (C.3) 

m=l m=l 

and 

4 4 . 

T^p = E^«1^' ^ E ^V;.(^«I^/3 + ^pD^)i^, . (C.4) 

i=l i=l 
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These operators represent the energy-momentum tensors for scalar and fermion fields, respec- 
tively. Under renormalization, they mix with the total energy-momentum tensor, which also 
includes the respective operator for the gluon fields and reads 

. 4 6 

T^^y = F^xK^ + ^uD^)^i + D^^rrMm + ■■■ 

i=l m=l 

^n, + T^, + T^^,. (C.5) 

Their mixing is governed by the anomalous dimension matrix for twist-two operators. The 
eigenoperators of the anomalous dimension matrix are 

Tr = T3 + + T^, Tii = -2T3 + T'^' + 2T't', Tm = -T^ + AT"^ - GT'^ . (C.6) 

The last two operators (unlike the former) have nonzero anomalous dimensions. After decom- 
posing the operators which appear in the OPE ( C.2 ) in terms of the above eigenvectors, i.e., 

12 11 

2 3 3 6' 

= + —Tii - —Tin , (C.7) 

3 15 21 35 ' ^ ' 



we finally get 

i I d^xe-^V'^(x) J-(0) = ^ (^v'-'^v"^ - ^r,-^ " r/^"^ + ^'''^) ^"/^ 

+ - , (C.8) 



15(Q 

from which one can read the coefficients of T^,y in the OPE of the current-current correlator. 
Although explicitly obtained here via a lowest-order calculation in perturbation theory, these 
coefficients turn out to be exactly as those (indirectly) computed at strong coupling, in Sect. 3. 
To see that, let us specialize (A^) to the high-energy regime, where 

q'^q^T^p ~ {q-fT^- ~ 2q^T^^ . (C.9) 

and then take the thermal expectation value by using the expression (A. 4) for the average energy 
momentum-tensor Q^^, = {T^u) in a strongly-coupled SYM plasma. We thus obtain 

^ [ d^xe--(J-(x)J^(0)) = (n-n^ - ^g-n^ - '-^g'^n^ + 



15(Q2)2 Q2 

TT^N^T^ f n-q \( n ■ q \ fifTV f (n 



6g2 Q2 '^/^y y-'^ Q'^^ ) io(Q2)2 qi 

which is in full agreement with (l3.20| ), as anticipated (recall that x = Q'^/2qT). Normally, the 
OPE coefficients at strong coupling are extracted by studying 3- and 4-point correlation func- 
tions. Our method in Sect. 3 is more straightforward (though limited to the energy momentum 
tensor) in that we do not have to compute higher point functions, but only use the known value 
of (T^u) at finite temperature. 
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D. High energy: the WKB approximation 



In this Appendix, we shah construct approximate solutions to the gravitational wave equations 
in the high energy regime at /c 3> ■ We shall thus confirm and extend the results found in 
Sect. 4, which, we recall, were valid only for u ^ 1. The complete solutions will be obtained by 
matching three different approximations, valid for different values of u: the two limited solutions 
valid for u ^ 1 and near n = 1, respectively, and the WKB solution valid in the intermediate 
region at < u <C 1. In this construction, the outgoing-wave condition will be imposed near 
the black hole horizon at u = 1, in conformity with the original prescription in Refs. pi] , |6^ . 

The general solutions valid for « <C 1 have been already constructed in Sect. 4. In the 
longitudinal sector, this is given by Eq. (|4.4| ), where the coefficient C2 is fixed by the boundary 
condition at n = 0, with the result shown in Eq. (|4.5| ); as for ci, this will be here obtained by 
matching onto the solution near n = 1, via the intermediate WKB solution. 

Consider now the solution near the horizon. For k ^ and u ~ 1, Eq. ( |2.14 ) simplifies to 



The solution which obeys the right, outgoing- wave, behaviour near u = 1 reads 

V;(n) =C3(l-n)^(i-*^). (D.2) 

(The second independent solution (1 — u) 2(^+*'^) must be rejected since it would describe a wave 
coming out from the horizon, i.e., a wave reflected by the black hole.) 

Furthermore, in the intermediate region <C it <C 1, with uq = 1/(4/;;^)^/^, the 'Schrodinger' 
wave equation reads 

iao(u) 



The WKB solution has the standard structure il^{u) = e^^om ^^|^/^| 



cjo(m) = / du'^/-V{u') = ±k du'- 
Jo Jo 1 



±'1 ( In - 2 arctan ] . (D.4) 
2 \ I- X u 



Hence the general solution in this intermediate region reads 



i^{u) = J ^-j=- [ciFk{u) + C5F_fc(u)], 

Fi^(u) = f ^^^ ^ ^ g-ifcarctanV^_ 

VI - VuJ 

We can now determine the unknown coefficients by matching the previous solutions in their 
common ranges of applicability. Comparing (p.2|) and (p.5|) near u = 1 gives 




C5 = 0, C3 = Ci\ll e^^''^2-*W4_ p_g) 
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Then a comparison of ( |4.4| ) and (D.5) in the region u <^ 1 but n S> — in this region 
ip{u) ~ y/ua{u) and ,^ ^ 1, so one can use the asymptotic expansions for the Bessel functions 
in Eq. (f4.4[) — gives, after simple calculations, 



ci 



C4 = A/-C2e 



-3j7r/4 



(D.7) 



As anticipated, we have recovered the simple relation (4^) between ci and C2 which implies that 
already the small-ti solution, Eq. ( [4.4D , is an outgoing wave, cf. Eq. ([4.7|). 

Turning now to the transverse sector, where the small-u solution was given in Eq. (^]^), we 
can similarly obtain the (outgoing-wave) solution near n = 1 as 



and the corresponding WKB solution as (compare to Eq. (p.5|) ) 



( ^~^ p-jfc arctan 



(D.9) 



(We have anticipated that C5 is set to zero after matching onto Eq. ( p.S] ).) The matching 
conditions then yield 



C4 = v3/7rci ( — I e 



I p 5j7r/12 



(D.IO) 



One finally gets the same result at small u as previously displayed in Eq. ([4.1C1|). 
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